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ABSTRACT 

We have studied a sample of 296 faint (> 0.5 mJy) radio sources selected from an area of 
the Tenth Cambridge (10C) survey at 15.7 GHz in the Lockman Hole. By matching this cata- 
logue to several lower frequency surveys (e.g. including a deep GMRT survey at 610 MHz, a 
WSRT survey at 1 .4 GHz, NVSS, FIRST and WENSS) we have investigated the radio spectral 
properties of the sources in this sample; all but 30 of the 10C sources are matched to one or 
more of these surveys. We have found a significant increase in the proportion of flat spectrum 
sources at flux densities below «1 mJy - the median spectral index between 15.7 GHz and 
610 MHz changes from 0.75 for flux densities greater than 1.5 mJy to 0.08 for flux densities 
less than 0.8 mJy. This suggests that a population of faint, flat spectrum sources is emerging 
at flux densities < 1 mJy. 

The spectral index distribution of this sample of sources selected at 15.7 GHz is compared 
to those of two samples selected at 1.4 GHz from FIRST and NVSS. We find that there is a 
significant flat spectrum population present in the 10C sample which is missing from the 
samples selected at 1.4 GHz. The 10C sample is compared to a sample of sources selected 
from the SKADS Simulated Sky by Wilman et al. and we find that this simulation fails to 
reproduce the observed spectral index distribution and significantly underpredicts the number 
of sources in the faintest flux density bin. It is likely that the observed faint, flat spectrum 
sources are a result of the cores of FRI sources becoming dominant at high frequencies. These 
results highlight the importance of studying this faint, high frequency population. 

Key words: galaxies: active - radio continuum: galaxies - galaxies: starburst 



1 INTRODUCTION 

The Tenth Cambridge Survey (10C; AMI Consortium: Franzen et 
al. 2011; AMI Consortium: Davies et al. 2011) at 15.7 GHz is the 
deepest high frequency (> 10 GHz) radio survey to date, complete 
to 1 mly in ten different fields covering a total of a 27 deg 2 . A 
further » 12 deg", contained within these fields, is complete to 
0.5 mJy. The 10C survey therefore enables us to study the faint 
source population at 15.7 GHz, a parameter space which has not 
been explored in any detail. Most studies of the sub-mjy popula- 
tion have focused on lower frequencies, due to the increased time 
required to survey a field to an equivalent depth at higher frequen- 
cies. 
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One area of particular interest has been the variation with flux 
density of the spectral index distribution for the samples selected 
at 15 - 20 GHz. For example, in the Ninth Cambridge survey (9C; 
Waldram et al. 2003) the proportion of sources with flat or rising 
spectra decreased as the flux density decreased. The median spec- 
tral index between 15.2 and 1.4 GHz (a| 5 4 2 ) changed from 0.23 for 
the highest flux density bin (S 15 ghz > 100 mly) to 0.79 for the low- 
est flux density bin (5.5 mJy > S 15 ghz > 25 mJy). (The convention 
5 oc v~ a , where 5 is flux density at frequency v, is used throughout 
this work). 

The Australia Telescope 20 GHz survey (AT20G; Massardi 
et al. 2011a) found a similar variation of spectral index with flux 
density for a sample with a flux density limit of 40 mly. The 10C 
survey (AMI Consortium: Davies et al. 2011) enables us to extend 
these studies to lower flux densities. Although the study by AMI 
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Consortium: Davies et al. (2011) contains a larger number of limits, 
it is clear that the fraction of flat spectrum sources increases again 
as flux density decreases further. 

There have been several attempts to model the high frequency 
radio sky (Wilman et al. 2008, Wilman et al. 2010, de Zotti et 
al. 2005, Tucci et al. 2011). These simulations have extrapolated 
from lower frequency surveys, generally on the assumption that the 
sources have power-law synchrotron spectra. These models are an 
increasingly poor fit to the observed sources counts below 10 mjy 
at 15 GHz. The number of sources is significantly underestimated, 
indicating that the properties of these sources are not well under- 
stood, largely due to the complexity and diversity of the high fre- 
quency spectra of individual sources. To understand the nature of 
the faint, high frequency population and constrain the models bet- 
ter, a multifrequency study is required. In this paper we describe 
just such a study, examining the radio properties of a sub-sample 
of 10C sources in the Lockman Hole, a region which has been ob- 
served over a wide range of wavelengths. 

The Lockman Hole is a region of the sky centred near 10 h 45 m , 
+58° (J2000 coordinates, which are used throughout this work) 
with exceptionally low HI column density (Lockman et al. 1986). 
The low infrared background (0.38 MJy sr~' at 100 |im; Lonsdale 
et al. 2003) in this area of the sky makes it ideal for infrared ob- 
servations. As a result, as part of the Spitzer Wide-area Infrared 
Extragalactic survey (SWIRE; Lonsdale et al. 2003) sensitive in- 
frared observations of « 14 deg 2 of the Lockman Hole area have 
been made. The availability of deep infrared observations in the 
Lockman Hole has triggered deep observing campaigns at optical, 
X-ray and radio wavelengths. 

The availability of data at such a wide range of frequen- 
cies makes the Lockman Hole a particularly good area for study. 
Here, the radio properties of sources detected in the 10C survey 
at 15.7 GHz in the Lockman Hole are investigated. The 10C data 
are combined with those from deep surveys at 610 MHz made with 
the Giant Meterwave Radio Telescope (GMRT; Garn et al. 2008, 
2010) and at 1.4 GHz made with the Westerbork Synthesis Ra- 
dio Telescope (WSRT; Guglielmino et al. 2012), along with other 
available data over a range of radio frequencies. The surveys used 
are described in more detail in Section 2 along with the process 
for matching them to the 10C catalogue. Section 3 describes how 
the parameters of the 10C sources are investigated using these sur- 
veys, including calculating the radio spectral indices and investig- 
ating the extent of the radio emission. The data are analysed and 
the properties of the 10C sources are discussed in Section 4. The 
10C source population is compared to samples selected at 1.4 GHz 
in Section 5 and to the Wilman et al. (2008, 2010) simulated model 
of the radio sky in Section 6. Optical identifications, redshift estim- 
ates and detailed discussions of the source types will be presented 
in a separate paper. 

Throughout this paper the term 'flat spectrum' refers to an ob- 
ject with spectral index a < 0.5 and 'steep spectrum' to an object 
with a > 0.5. 



2 SAMPLE SELECTION 
2.1 Surveys used 

The 10C radio survey at 15.7 GHz was made with the Arcminute 
Microkelvin Imager (AMI; Zwart et al. 2008) with a beam size of 
30 arcsec. It covers «27 deg 2 complete to 1 mJy and «12 deg 2 
complete to 0.5 mJy across ten different fields; a full description 



of the 10C survey and the source catalogue can be found in AMI 
Consortium: Franzen et al. (201 1) and AMI Consortium: Davies et 
al. (201 1). Two of the fields are in the Lockman Hole, covering an 
area of 4.64 deg 2 (see Fig. 1) and detecting a total of 299 sources. 

To investigate the properties of these 10C sources we have 
matched the 10C catalogue to other lower-frequency, but usually 
higher-resolution, radio catalogues as detailed below. This not only 
enables us to determine the radio spectral properties of the sources 
but also allows us to investigate the extent and structure of the ra- 
dio sources in more detail. The greater positional accuracy of these 
higher-resolution catalogues is also vital for finding the counter- 
parts of the 10C sources at optical and infrared wavelengths; the 
analysis of the optical and infrared properties of the 10C sources 
will be presented in a later paper. The other radio surveys of the 
Lockman Hole used in this work are listed in Table 1 and briefly 
described below. 

A series of deep observations at 610 MHz made with the 
GMRT covers the whole 10C Lockman Hole area except for a small 
corner of the field containing five 10C sources (Fig. 1). The GMRT 
image has an rms noise of «60 uJy per beam in the central area (see 
Garn et al. 2008, 2010, for details of the data reduction and source 
extraction). This deep image and the catalogue derived from it are 
used here. 

A deep survey at 1.4 GHz carried out with the WSRT overlaps 
a large portion of the 10C survey (Fig. 1) in the Lockman Hole 
(Guglielmino et al. 2012). The rms noise in the centre of the map 
is a; 1 1 |ijy; this map and the associated source catalogue are used 
in this study. 

We also make use of several other catalogues and surveys cov- 
ering the Lockman Hole region. The Faint Images of the Radio Sky 
at Twenty cm (FIRST; White et al. 1997), NRAO VLA Sky Survey 
(NVSS; Condon et al. 1998) and Westerbork Northern Sky Sur- 
vey (WENSS; Rengelink et al. 1997) cover the whole area. There 
are also several deep observations made with the VLA within the 
Lockman Hole area: Owen et al. (2009) at 324 MHz (OMK2009) 
and Biggs & Ivison (2006) and Owen & Morrison (2008) at 1.4 
GHz (BI2006 and OM2008). The locations of these surveys are 
shown in Fig. 1 and are summarised in Table 1. 

Sixteen 10C sources in the Lockman Hole are labelled in the 
10C catalogue as part of a 'group' (see AMI Consortium: Davies 
et al. 201 1 for details) - in this case consisting of eight pairs. Con- 
tour maps of these eight pairs were examined by eye, along with 
images of their counterparts in FIRST and/or GMRT. For three of 
these pairs, there is evidence of structure connecting the two com- 
ponents, so it was decided to combine the two components into one 
source (an example of one such source is shown in Fig. 2). The po- 
sitions listed in the catalogue for these three pairs is the point mid- 
way between the two components. For the remaining five pairs the 
components were left as separate sources. This leaves a total of 296 
sources in the 10C sample. 

These 296 sources, selected at 15.7 GHz, are the subject of this 
paper. A subsample of 89 sources, with flux densities greater than 
0.5 mJy, selected from the region complete to 0.5 mJy is defined 
as 'Sample A'. A second subsample of 118 sources, 'Sample B', 
which has some overlap with Sample A, forms a sample complete 
to 1 mJy. Sixty- two 10C sources, 43 of which form a sample com- 
plete to 0.5 mJy at 15.7 GHz, are in the deep areas observed with 
the VLA at 1.4 GHz by Biggs & Ivison and Owen & Morrison 
(B 12006 and OM2008). All 62 sources are detected at 1.4 GHz 
so spectral information is available for all of them. They therefore 
form a particularly useful subsample, defined as Sample C. Table 2 
contains a summary of the different subsamples used in this paper. 
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Figure 1. The deep radio surveys in the Lockman Hole region. The two lar- 
ger (red) squares show the shallow region of the IOC survey fields, complete 
to 1 mJy. The two smaller squares contained within the larger squares are 
the deep IOC regions, complete to 0.5 mJy. The pentagon (dark blue) shows 
the GMRT survey area, the dashed square (green) shows the OM2008 sur- 
vey area, the smaller circle (pale blue) shows the BI2006 survey area and 
the larger circle (pink) shows the WSRT survey area. FIRST, NVSS and 
WENSS are not shown as they cover the whole region. See Table 1 for 
details of the different surveys shown. 



2.2 Matching the radio catalogues 

2.2.7 Choosing a match radius 

The Topcat software package 1 was used to match the catalogues. 
The match radius was chosen so as to maximise the number of real 
associations and avoid false matches. The error in the 10C posi- 
tions is a; 6 arcsec, which is larger than the errors in the other cata- 
logues used here (for example the error in the GMRT positions is 
< 1 arcsec), so the errors in the 10C positions tend to dominate. 
Here we describe the process by which we chose a suitable match 
radius and assessed the probability of genuine and random matches 
for the 10C and GMRT catalogues. To create a random distribution 
of sources the 10C sources were shifted by 5 arcmin in declina- 
tion. Both this shifted catalogue and the true 10C catalogue were 
then matched to the GMRT catalogue and the angular separation 
between a shifted or true 10C source and its nearest GMRT source 
was recorded. Once a shifted or true 10C source was matched to 
a source in the GMRT catalogue no further matches were sought 
for that source. The matches were then binned according to their 
separation for both the true and shifted sources. The resulting dis- 
tribution is shown in Fig. 3. It can be seen that beyond 15 arcsec the 
numbers of random matches become comparable for both the true 
and shifted catalogues; thus 15 arcsec was chosen as the match ra- 
dius. The number of matches with the shifted distribution within a 
separation of 15 arcsec is 4 (out of 296 sources), meaning the prob- 
ability of a random match within 15 arcsec is a; 1.5 percent of the 
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Figure 2. 15.7-GHz image of an example extended source listed as two 
components in the 10C catalogue. The contours are plotted at (±2 V2", n = 
0, 1...7) X 0.136 mJy. The crosses mark the positions of the two sources in 
the 10C catalogue. There are three such pairs of sources in the 10C Lock- 
man Hole fields; in each case the flux densities of the separate components 
are combined. 
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Figure 3. The number of matches with a given separation between the 10C 
and GMRT catalogues for the true and shifted distribution of 10C sources. 
Beyond 15 arcsec the number of true and shifted matches becomes compar- 
able so 15 arcsec is chosen as the match radius. 



total number of sources being matched. 15 arcsec corresponds to 
2.5cr, where a is the typical error in the 10C source positions, sug- 
gesting that the matching is limited by the error in the 10C source 
positions. 

This process was repeated for the other catalogues and it was 
decided that 15 arcsec was a suitable match radius to use when 
matching each of the catalogues to the 10C catalogue. The number 
of random matches within 15 arcsec for each catalogue is shown in 
Table 3. 
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Table 1. Radio catalogues in the Lockman Hole 



Catalogue 


Reference(s) 


Epoch of observation 


Frequency 


Beam size 


rms noise 


Area covered 








/GHz 


/arcsec 


/mJy 


/deg 2 


1 0r* shallow 4 


AMI Consortium: Franzen et al. (2011) 
AMI Consortium: Davies et al. (2011) 


AiiiT 7008 - Time 7010 


15.7 


30 


0.1 


4 (\A- (f\ppr\ nrpa*; itipIiiHpH^ 


IOC - deep I 


AMI Consortium: Franzen et al. (201 1) 
AMI Consortium: Davies et al. (2011) 


Aug 2008 - June 2010 


15.7 


30 


0.05 


1.73 


GMRT i 


Gam et al. (2008) 
Garnet al. (2010) 


Jul 2004 - Oct 2006 


0.610 


6x5 


0.06 


13 


WSRT 


' Guglielmino et al. (2012) 


Dec 2006 - Jun 2007 


1.4 


11x9 


0.011 


6.6 


OM2008 


Owen & Morrison (2008) 


Dec 2001 -Jan 2004 


1.4 


1.6 


0.0027 


0.011 


OMK2009 


Owen et al. (2009) 


Feb 2006 - Jan 2007 


0.324 


6 


0.07 


3.14 


BI2006 


Biggs & Ivison (2006) 


Jan 2001 - Mar 2002 


1.4 


1.3 


0.0046 


0.089 


FIRST 


White et al. (1997) 


1997-2002 


1.4 


5 


0.15 


Whole area 


NVSS 


Condon et al. (1998) 


1997 


1.4 


45 


0.45 


Whole area 


WENSS 


Rengelink et al. (1997) 


1991 - 1996 


0.325 


54 


3.6 


Whole area 



Table 2. The different subsamples in the Lockman Hole used in this paper 



Sample 


Description 


Number of sources 


All sources 


All sources detected at 15.7 GHz. Includes sources below 0.5 mJy. 


296 


Sample A 


Complete to 0.5 mJy at 15.7 GHz 


89 


Sample B 


Complete to 1.0 mJy at 15.7 GHz 


118 


Sample C 


All sources in the deep regions surveyed at 1 .4 GHz with the VLA by Owen & Morrison 
(2008) and Biggs & Ivison (2006). This sample contains sources below the completeness 
limit at 15.7 GHz. 


62 



Table 3. A summary of the matching between the 10C catalogue and other 
radio catalogues in the Lockman Hole. Column (2) shows the number of 
matches to 10C sources within 15 arcsec and column (3) shows the number 
of 10C sources in the survey area. Column (4) shows the number of random 
matches within 15 arcsec, i.e. the number of matches to the shifted 10C 
sources. 



Catalogue 


Number of 


Number of 


Number of 




matches 


10C sources 


random matches 


(1) 


(2) 


(3) 


(4) 


GMRT 


205 


291 


8 


WSRT 


160 


182 


10 


FIRST 


196 


296 





NVSS 


166 


296 





WENSS 


86 


296 


6 


OM2008 


27 


27 


2 


OMK2009 


116 


156 


2 


BI2006 


35 


35 


3 



Table 4. A summary of the number of 10C sources with multi-frequency 
matches. 



Number of matches 


Number of sources 


10C only 


30 


One match only 


46 


Two matches 


36 


Three matches 


27 


Four matches 


55 


Five matches 


78 


Six matches 


23 


Seven matches 


1 



3 DERIVING THE SOURCE PARAMETERS 



2.2.2 Results of the matching 

A summary of the matching between the 10C sources and the other 
radio catalogues is shown in Table 3. In total, 266 out of the 296 
10C sources are detected at at least one other wavelength. Table 4 
summarises the number of sources with matches in more than one 
other catalogue. 

For the 30 10C sources with no other detections, upper limits 
on the flux density at 610 MHz from GMRT and at 1.4 GHz from 
WSRT and FIRST can be used, as described in Section 3.1. All of 
the 10C sources with no other detections have a low flux density 
(see Fig. 4), with the majority below 1 mJy and all below 2 mJy. 



3.1 Flux density values 

3.1.1 GMRT 

The large difference in resolution between the 10C (beam size 30 
arcsec) and GMRT (beam size 6 arcsec) observations means that 
care must be taken when comparing flux densities. As illustrated 
in Fig. 5, many 10C sources are extended or resolved into mul- 
tiple components at 610 MHz. The components of such sources 
will be listed as separate entries in the GMRT catalogue. In addi- 
tion, for some of the sources in the GMRT catalogue there may be 
extended low brightness structure too faint to be seen in the high 
resolution images. In order to reduce these problems, contour plots 
of all the GMRT sources which matched to a source in the 10C 
catalogue (using a 15 arcsec match radius as described in Section 
2.2.1) were examined by eye. For the 50 sources which appeared 
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Figure 4. The distribution of 10C sources with and without associations in 
the other radio catalogues as a function of flux density. This diagram con- 
tains all 296 sources detected, some of which are below the completeness 
levels of the 10C survey. All of the 10C sources with no matches have a flux 
density below 2 mJy. 

multiple or extended, the GMRT images were convolved in AIPS 2 
to a 30-arcsec gaussian to create an image of comparable resolu- 
tion to the 10C data (see Fig. 5). The total 610-MHz flux density of 
each source (which had already been matched to a 10C source from 
the higher resolution image) was then estimated from the smoothed 
GMRT image. The integrated 610-MHz flux densities were estim- 
ated by fitting a Gaussian using the AIPS task JMFIT. In the few 
cases where JMFIT did not converge (due to the presence of an- 
other bright source in the subimage) the integrated flux density was 
found by hand using the AIPS task TVSTAT. 

For the 84 10C sources for which a counterpart was not 
present in the full resolution GMRT catalogue the GMRT images 
were used to place an upper limit on the source flux density at 
610 MHz. A sub image of 2.5 arcmin was extracted and smoothed 
(as described above) to check for any large scale structure which 
might have been resolved out in the original image. In eight cases 
a source could be seen in the smoothed GMRT images with a peak 
within 15 arcsec of the 10C source position; for these sources a 
value for the integrated flux density was obtained manually using 
TVSTAT. For the remaining sources, the upper limit on the flux 
density was taken to be three times the noise in this smoothed im- 
age. It was not possible to get any information about the flux dens- 
ity at 610 MHz for the five 10C sources outside the GMRT area 
so these sources are excluded from all discussions relating to 610- 
MHz data. 

3.1.2 NVSS 

NVSS has a resolution comparable to the 10C survey so the flux 
densities can be directly compared without any problems caused 
by resolution differences. Because of this, it is useful to be able 
to calcuate spectral indices using NVSS and 10C values only, so 
upper limits were found for all sources not matched in the NVSS 
catalogue despite the fact that some of these sources have a coun- 
terpart in one of the deeper 1.4-GHz surveys used. NVSS images 

2 Astronomical Image Processing System, see http://www.aips.nrao.edu/. 



of all IOC sources without a match in the NVSS catalogue were 
examined to see if there was a source present which was below the 
catalogue limit and for the few sources where this was the case, 
the flux density of the NVSS source was found manually using the 
AIPS task TVSTAT. For the remaining unmatched sources, an up- 
per limit of three times the local noise in the NVSS image was 
placed on the flux density. 

3.1.3 WSRT and FIRST 

As the WSRT synthesised beam (12 arcsec) is smaller than the 
IOC synthesised beam (30 arcsec), several sources were resolved 
into multiple components in the WSRT catalogue but appeared as 
a single component in the 10C catalogue. In order to reduce the ef- 
fects of resolution when comparing the flux densities, sub-images 
of the WSRT image at the position of each 10C source were ex- 
amined by eye; for those sources where there were multiple WSRT 
components associated with one 10C source, the flux densities of 
the WSRT components were added together. (In fact, in all cases 
except one, the WSRT sources in question were listed as a multi- 
component source in the WSRT catalogue and the total flux density 
from the catalogue was used). This process was also carried out for 
FIRST sources. 

For those 10C sources in the WSRT survey area without a 
match to the WSRT catalogue, the WSRT image was examined to 
see if there was a source present which was below the catalogue 
limit. This was not the case for any of the unmatched 10C sources 
in the WSRT area. The upper limit of flux density was taken to be 
three times the local noise in the WSRT image. For the 28 sources 
which were unmatched at 1.4 GHz and outside the WSRT area, an 
upper limit of 1 mjy was placed on the flux density because this is 
the FIRST completeness limit and these sources are all within the 
FIRST survey area but not detected. 

3.1.4 OM2008 and BI2006 

The two deep 1.4 GHz surveys made with the VLA (OM2008 and 
BI2006) have very small synthesised beams (« 1.5 arcsec) com- 
pared to the 10C survey. As images are not available for these 
surveys, the flux density values from the catalogues were used. 
OM2008 do account for the fact that some of the sources may be 
extended when calculating flux densities. They convolve the full 
resolution images to an effective beam size of 3, 6 and 12 arc- 
sec and compare the flux densities derived from the four images. 
BI2006 do not attempt to account for extended sources when cal- 
culating the flux densities, however there is only one 10C sources 
which has a counterpart in BI2006 and not in any other 1 .4-GHz 
catalogue so this will not significantly affect these results. 

3.2 The effects of variability 

The different surveys used in this paper were not carried out simul- 
taneously so it is important to consider the possible effects of vari- 
ability on the observed spectral index distributions. The epoch of 
observation for each survey is shown in Table 1. The time interval 
between the 15-GHz and the 610-MHz observations is in the region 
of 4 - 6 years and between the 15-GHz and 1.4-GHz observations 
5-10 years. 

Whilst there is currently no data on the variability of sources 
at 15 GHz at the low flux density end of our sample, there have 
been some systematic studies at higher flux densities. Bolton et al. 
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Figure 5. The difference in resolution between AMI Large Array (used for the IOC survey) and the GMRT. Figs, (a) and (d) each show a source imaged with 
AMI, Figs, (b) and (e) show the same sources imaged with GMRT. In both cases several sources in the GMRT catalogue correspond to one source in the IOC 
catalogue, so care must be taken when comparing the catalogues. In Figs, (c) and (f) the GMRT images have been smoothed to create an image with a resolution 
comparable to the IOC images. This allows for a more direct comparison between the two surveys. The contours are drawn at (±2 V2", n = 0, 1...7) X x mJy 
where x = 0.074 for (a), 0.1 for (b), 0.3 for (c), 0.086 for (d), 0.16 for (e) and 0.3 for (f). 



(2006) studied the 15-GHz variability of 51 9C sources with flux 
densities > 25 mJy over a 3-year period. They found that while 
there was no evidence for variability (above the ~6 percent flux 
calibration uncertainties) in steep spectrum sources, half of the flat 
spectrum objects were variable. In total, 29 percent of the sources 
studied were found to vary. Sadler et al. (2006) observed 173 20- 
GHz sources with flux densities > 100 mJy over a 2-year period 
and found that 42 percent varied by more than 10 percent. How- 
ever, they found no correlation between variability and radio spec- 
tral index. More recently, Bonavera et al. (2011) investigated the 
variability of 159 sources with 20-GHz flux densities > 200 mJy 
and found the variability to be slightly larger than that found by 
Sadler et al., with an r.m.s. amplitude of 38 percent at 20 GHz on 
a timescale of a few years. Bonavera et al. note that there is some 
indication that the variability decreases as flux density decreases. 
Massardi et al. (201 lb) studied the variability of a brighter sample 
(52oghz > 500 mJy) at 20 GHz and found similar levels of vari- 
ability to Sadler et al. At higher flux densities still Franzen et al. 
(2009) have looked at variability at 16 GHz in a complete sample 
of 97 sources with flux densities > 1 Jy over timescales of about 
1.5 years. They found that 15 percent of the sources vary by more 
than 20 percent; however, in contrast to the results of Sadler et al. 
but in agreement with those of Bolton et al., the spectra of the vari- 
able sources are flatter than those of the non-variable ones. 

The variability properties of the population studied here are 



not known. However, on the assumption that the faint sources in the 
10C sample exhibit the same sort of flux densiy variations as shown 
in the higher flux density samples discussed above, a significant 
fraction of them are likely to have varied over the period between 
the observations. Thus the spectral indices of individual objects 
may be unreliable. However, given that the sources are probably 
equally likely to increase or decrease in flux density this should not 
have a major effect on the overall spectral index distribution. 

3.3 Spectral indices 

To investigate the spectral properties of the source sample in 
a quantitative way, the spectral index was calculated between 
15.7 GHz and 1.4 GHz (aj 5 4 7 ) and 15.7 GHz and 610 MHz (a' 5 J) 
for each source. For a| 5 4 7 all 296 sources are studied, for a-^ 7 the 
5 sources outside the GMRT area are excluded as there is no 610- 
MHz flux density information available. For the sources with no 
match in GMRT, a limiting spectral index was calculated from the 
upper limit placed on the flux density from the GMRT image, as 
described in Section 3.1.1. 

The distributions of a] 5 4 7 were investigated in two ways, using 
slightly different procedures, to check the effects of resolution on 
the data. The first, a\\ 7 M, makes use of all of the 1.4-GHz data 
available; for the sources where there is more than one 1.4-GHz 
flux density, flux densities are chosen according to resolution in 
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Figure 6. Compactness ratio calculated using GMRT data, Rqmrt, against 
the signal to noise ratio from the GMRT catalogue. The vertical line shows 
the divide at peak flux density / noise = 20 and the horizontal lines at R = 
1.35 and 2.00 indicate the cutoffs in R used when classifying the extended 
sources, see text for details. 



Figure 7. Value of the compactness ratio, R, calculated using FIRST and 
GMRT data. The horizontal and vertical lines indicate the values used when 
identifying the extended sources; see text for details. Sources classified as 
extended and those not classified are shown separately. Sources with a sig- 
nal to noise ratio < 20 in the GMRT data are indicated with a green circle. 



the following order of preference: NVSS, WSRT, OM2008, FIRST, 
BI2006; FIRST and BI2006 are last as they are the most likely to 
resolve out some of the flux of the 10C sources because of their 
small beam sizes (5 and 1.5 arcsec respectively). OM2008 also 
has a small beam (1.6 arcsec) but the sources have been convolved 
with gaussians of varying radius to try and overcome the resolution 
problem. NVSS and WSRT have larger beam sizes (30 and 12 arc- 
sec respectively) which are more comparable to the 10C beam. For 
the sources with no match in any of the 1 .4-GHz catalogues which 
are in the WSRT survey area, the upper limit from the WSRT im- 
age (as described in Section 3.1.3), is used to calculate a limiting 
spectral index. For the remaining sources, an upper limit of 1 mjy 
from the FIRST survey is used. 

The second value, a\ 5 ^N, only uses values from the NVSS 
catalogue as this has a resolution comparable to the 10C survey. 
For the 10C sources which do not appear in NVSS, the limit is de- 
rived from the NVSS image (see section 3.1.2). These values of the 
spectral index contain a larger number of upper limits but provide 
a useful comparison when considering the effects of resolution. 

Possible effects of variability on the spectral indices are dis- 
cussed in Section 3.2. 

3.4 Extent of the radio emission 

To determine in a quantitative manner whether a source is extended 
or not, a compactness ratio R is often used. This is usually taken as 
the ratio of the integrated or total flux density of a source 5 ,„, and 
its peak flux density on a map 5,,^ i.e. R = S in ,/S pea k- For the 
10C sample considered here we do not use the 10C observations to 
find R because the large beam size of « 30 arcsec means that the 
majority of the sources are unresolved. Instead, we use the matched 
data from the lower frequency catalogues, which have smaller beam 
sizes and can therefore provide more information about the angular 
size of a source. Four values of R are calculated, using data from 
the GMRT, FIRST and OMK2009 (324 MHz VLA observations) 
catalogues, which all have beam sizes ss 6 arcsec, and the WSRT 
catalogue which has a beam size of 12 arcsec. The deep 1.4 GHz 
VLA observations (OM2008 and BI2006) are not used in this ana- 



lysis because their beam sizes are considerably smaller (as 1 arcsec) 
so that the data are not comparable with those from the other cata- 
logues used here. 

To take account of the effects of noise, several other studies 
(e.g Bondi et al. 2007, Prandoni et al. 2006) plot R against the sig- 
nal to noise ratio and fit a lower envelope to the data; reflecting 
this about R = 1 gives a curve which provides the cutoff between 
the extended and compact sources. We have not used this method 
for the sample of sources in this paper due both to the relatively 
small number of sources in the sample and to the range of differ- 
ent surveys used. Instead, each source was examined by eye and the 
criteria described below were decided upon to identify the extended 
sources. 

At lower signal to noise ratios the errors in the R values be- 
come larger, due to the larger errors in the integrated and peak flux 
density values. This is evident in Fig. 6 for the GMRT data, where 
the number of sources with values of R < 1 , which must be due to 
errors in the flux density values, increases at lower signal to noise 
ratios. Therefore different criteria are used to identify the extended 
sources at high and low signal to noise. Sources were classified as 
extended if they had a signal to noise ratio greater than 20 in the 
GMRT data and values of R > 1.35 in at least two of the datasets 
used, or a value of R > 2 in at least one of the datasets. Sources 
with a signal to noise ratio less than 20 in the GMRT data were 
classified as extended if the value of R was greater than 2 in any 
of the datasets. A comparison of two of the R values used is shown 
in Fig. 7. For most sources the two values are similar and some of 
the largest variations are for those sources with low signal to noise 
values, which is expected because the errors in the integrated flux 
densities are larger for the fainter sources. 

We are confident that all the sources which fulfil these criteria 
are extended; however, as we have erred on the side of caution, 
there will be some sources which have not been classified as exten- 
ded but do in fact have extended emission on the scale investigated 
here. In particular, sources which only have a value of R from the 
WSRT catalogue are not classified as extended as the WSRT beam 
is larger than that of the other catalogues used here. For this reason, 
those sources which are not classified as extended are placed in 
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Figure 9. The spectral index ff]^ 7 M against 15.7-GHz flux density, show- 
ing which catalogue the 1 .4-GHz flux density values come from. 
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Figure 8. Example spectra. 10CJ 104724+573703 is an example of a steep 
spectrum object, 10CJ104710+582821 is a rising spectrum object and 
10CJ105535+574636 is an example of a flat spectrum object. The blue 
cross at 324 MHz (if present) is from WENSS, the red plus at 324 MHz 
is from OMK2009 and the value at 610 MHz is from GMRT. The blue 
cross at 1 .4 GHz is from FIRST, the green circle is from WSRT and the red 
plus is from NVSS. The value at 15.7 GHz is from 10C. Error bars are not 
shown when they are smaller than the symbol plotted. 



1.5 

1 - 



x No limits 
v 1.4-GHz limits 

610-MHz limits 

1 .4-GHz and 61 0-MHz limits 



0.5- 


-0.5 

-1 - 
-1.5- 



'SI 



-1.5 



-0.5 



a 15 ' 7 

"fj.61 



0.5 



1.5 
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between steep and flat spectrum sources; sources with a > 0.5 are clas- 
sified as steep spectrum while sources with a < 0.5 are classified as flat 
spectrum. The dotted line illustrates al: 



,15.7 
0.61 



the 'unclassified' bin. The 36 sources which do not have a coun- 
terpart in any of the four catalogues used here are classified as 'no 
information'. The criteria used to select the extended sources here 
is roughly equivelent to selecting sources with angular sizes larger 
than « 6 arcsec. 

A catalogue containing the flux density, spectral index and R 
values will be available online. 



4 SAMPLE ANALYSIS 

4.1 Radio spectral properties 

Some example radio spectra are shown in Fig. 8, demonstrating the 
different spectral types observed. Spectral indices are calculated for 
all sources, as described in Section 3.3. The spectral index a\ s ^M 
against 15.7-GHz flux density for all sources is shown in Fig. 9. It is 
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Figure 10. The median spectral index and percentage of sources with a > 0.5 as a function of 15.7-GHz flux density for all 10C sources. The values of a [ [ 4 M 
are calculated using the best 1.4-GHz flux density values available (order of preference = NVSS, WSRT, OM2008, FIRST, BI2006) while the values of a\ 5 ?N 
are calculated using 1.4-GHz flux density values from NVSS only. The points are plotted at the median flux density in each bin but are offset for clarity. 
Limiting spectral indices are included in the median using survival analysis (see text). The error bars in the median spectral index plot show the interquartile 
range divided by to give an indication of the errors in the medians. Note the large change in spectral index between 1 and 2 mjy. 
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Figure 11. The median spectral index and percentage of sources with a > 0.5 as a function of 15.7-GHz flux density for sources in the complete samples 
A and B. The values of a]^ 7 M are calculated using the best 1.4-GHz flux density values available (order of preference = NVSS, WSRT, OM2008, FIRST, 
BI2006). The points are plotted at the median flux density in each bin but are offset for clarity. The sources in the lowest two flux density bins are from Sample 
A, which is complete above 0.5 mjy and the remaining three flux density bins contain sources from sample B, complete above 1 mjy. Limiting spectral indices 
are included in the median using survival analysis (see text). The error bars in the median spectral index plot show the interquartile range divided by ^fN to 
give an indication of the errors in the medians. 



clear from this plot that there is a greater proportion of flat spectrum 
sources at lower flux densities. This trend is further investigated by 
calculating the median spectral indices cr^ 7 and a\ 5 4 J M in three 
different flux density bins containing equal numbers of sources and 
the results are shown in Tables 5 and 6. Upper limits are included by 
using the ASURV Rev 1.2 package which implements the survival 
analysis methods presented in Feigelson & Nelson (1985). The dis- 
tributions of both a-Q 5 J and a\ 5 4 M as a function of 15.7-GHz flux 
density show a sharp change at flux densities between 1 and 2 mjy. 
The spectral index distributions for these three flux density bins for 
both or!" and a[ 5 4 7 M are shown in Fig. 13 and are very similar. 
There is a distinct peak at a ~ 0.7 in the highest flux density bin 
(Shigh). As flux density decreases the peak broadens as the contribu- 
tion from flat spectrum sources becomes much more significant. In 



the lowest flux density bin (Si ow ) the sources display a wide range 
of spectral index values, with a broad peak at a ~ 0.3. The sources 
shown in white are upper limits and could only move to the left in 
these plots, making the Si ow distribution even more different from 
the Shi e h distribution. 

This is illustrated in Fig. 10 which shows the median aJ 5 J, 
cr] 5 4 7 M and a\ 5 4 N and the percentage of sources with a > 0.5 in 
narrower 15.7-GHz flux density bins - there is a far higher propor- 
tion of flat spectrum sources at lower flux densities. The values for 
as] FN (using NVSS flux densities only) are very similar to those 
for a\ 5 4 7 M, except in the lowest flux density bin which is dominated 
by upper limits in the NVSS-only case. This implies that resolution 
differences between 10C and the other surveys is not having a ma- 
jor effect on the derived spectral indices. The same plots for the 
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Table 5. The aj" results for three 15.7-GHz flux density bins. 



Bin 


15.7-GHz flux density 


Number 


Number of 


Median a 


Mean a 


% a > 0.5 


name 


range / mJy 


of sources 


upper limits 








S low 


0.300 < S < 0.755 


99 


45 


0.08 


0.25 ± 0.04 


25 ±4 


^med 


0.755 < 5 < 1.492 


97 


29 


0.36 


0.31 ±0.05 


40 ±5 


S high 


1.492 < 5 < 45.700 


95 


4 


0.75 


0.57 ± 0.05 


67 ±5 



Table 6. The a\ 5 ^ M results for three 15.7-GHz flux density bins. 



Bin 


15.7-GHz flux density 


Number 


Number of 


Median a 


Mean a 


% a > 0.5 


name 


range / mJy 


of sources 


upper limits 








Slow 


0.300 < S < 0.755 


99 


27 


0.10 


0.11 ±0.07 


29 ±4 


^med 


0.755 <S < 1.492 


99 


13 


0.43 


0.30 ± 0.07 


46 ±5 


^high 


1.492 <S < 45.700 


98 


2 


0.79 


0.66 ± 0.05 


76 ±4 



complete samples A and B are shown in Fig. 11. There are fewer 
sources in this sample so the uncertainties are larger but the general 
trend towards decreasing spectral indices below as 2 mJy remains 
the same. 

A plot of a] 5 4 7 M against ai 5 ;] and for all sources is shown 
in Fig. 12. There is a good correlation between the two values of 
spectral index. The majority of the points in the bottom left corner 
which deviate from the correlation are upper limits and therefore 
could move closer to the one-to-one correlation line shown. Of the 
points which are not upper limits, slightly more lie above the one- 
to-one correlation line than below it, indicating a slight spectral 
steepening at higher frequencies. 

Possible effects of variability are discussed in Section 3.2. The 
variability properties of the population studied here are not known, 
but it seems likely that, if anything, variability would increase the 
proportion of steep spectrum sources rather than producing the de- 
crease observed. It is possible, for example, that variability in some 
of the genuinely flat spectrum sources over the epochs covered by 
the surveys used here could make them appear to have steep spec- 
tra. On the other hand, the genuinely steep spectrum sources in our 
faint sample are unlikely to be variable (unless their properties dif- 
fer markedly from the steep spectrum sources at higher flux densit- 
ies) so their spectral indices will not be affected at all. 
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Figure 14. The median spectral index as a function of flux density from the 
study of 9C sources by Waldram et al. (2010) and the 10C sources studied in 
this paper. The points are plotted at the bin midpoint, except for the highest 
flux density bin (S > 100 mJy) where the point is plotted at the bin lower 
limit. 



4.2 Comparison with other spectral index studies 

The variation in spectral index with flux density at higher (~ 10 mJy 
- 1 Jy) flux densities has been investigated by Waldram et al. (2010) 
and Massardi et al. (201 la). Waldram et al. and Massardi et al. used 
samples selected at 15 and 20 GHz respectively and, by matching 
them to catalogues at 1 .4 GHz, found that the median spectral index 
becomes rapidly larger with decreasing flux density; for example 
Waldram et al. found that the median a\\ changed from around 
at 1 Jy to 0.8 at 10 mJy. Our study of 10C sources shows that as 
the flux density decreases further the median spectral index drops 
again, with the median a[ 5 4 decreasing from around 0.8 for flux 
densities > 1.5 mJy to around 0.1 at 0.5 mJy. As discussed in Sec- 
tion 4.1 the change occurs fairly abruptly at ss lmjy, indicating that 
the nature of the 15-GHz source population is changing at this flux 
density. The relationship between spectral index and flux density is 
summarised in Fig. 14 which shows the Waldram et al. results from 
the Ninth Cambridge survey (9C) along with the results from this 
study. 

We have shown here that the spectral index distribution 



changes dramatically below ~1 mJy for a sample of sources se- 
lected at 15.7 GHz. A similar trend has been found by other studies 
at slightly lower frequencies. For example, Prandoni et al. (2006) 
studied a complete sample of 1 3 1 sources detected at 5 and 1 .4 GHz 
with a comparable flux density range to the sample studied here. 
They found that the median spectral index between 5 and 1.4 GHz 
changed from 0.56 ± 0.06 for sources with 5 5G hz > 4 mJy to 
0.24 ± 0.06 for sources with Ssghz < 4 mJy. 



4.3 Extent of the radio emission 

The sources were split into three groups (extended, unclassified 
and no information) as described in Section 3.4. A summary of 
the properties of these three groups of sources in terms of spectral 
index and flux density is given in Table 7 and the flux density distri- 
bution for the three groups is shown in Fig. 15. A larger proportion 
of the brighter sources are extended, with extended sources mak- 
ing up 38 percent of the total number of sources with S > 1 mJy, 
compared to 19 percent of sources with S < 1 mJy. This could be 
a result of the more stringent criteria for classifying sources with 
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Figure 13. The spectral index distribution for three different flux density bins are shown in the top three panels. Si ow : 0.300 < S15.7GHZ =S 0.755 mJy, 
Smed : 0.755 < S]5.7GHz < 1.492 mJy, Shigh : 1.492 < S15.7GHZ < 47 mJy. The spectral index distribution for the whole 10C sample is shown in the bottom 
panel. o^ 5 J is calculated using the GMRT values and limits, accounting for the resolution difference as described in Section 3.1, and using the integrated flux 
density from the 10C catalogue. a\ s 4 7 was calculated using the best 1.4 GHz flux density available (order of preference = NVSS, WSRT, OM2008, FIRST, 
BI2006). In both cases upper limits on a are included for sources where no low frequency flux density data is available; these are shown in white while values 
are plotted in black. 
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Table 7. A summary of the properties of sources classified as extended, 
those not classified and those with no information. The value of spectral 
index used here is al jj so the five sources which are outside the GMRT 
survey area are not included. The numbers in brackets refer to the percent- 
age of sources in each classification. 







Number of sources with 






Bin 


Total 


S 15 7 > 1 mjy 


S 15 7 < 1 mjy 


a > 0.5 


a < 0.5 


Extended 


85 


55 (65) 


30 (35) 


70 (82) 


15 (18) 


Unclassified 


170 


80 (47) 


90 (53) 


56 (33) 


1 14 (67) 


No info 


36 


5(14) 


31 (86) 


2(6) 


34 (94) 




0.5 1 
log(1 5.7-GHz flux density / mjy) 

Figure 15. The 15.7-GHz flux density distribution for the extended sources, 
those which are unclassified and those with no information. The three his- 
tograms are overlaid. 



low signal to noise as extended as well as the fact that extended 
low surface brightness emission is more likely to be missed at low 
signal to noise levels. As expected, the majority of the sources with 
no information are faint, with S < 1 mjy. 

The spectral index distributions for the extended sources and 
those not classified are, as anticipated, significantly different, as 
shown in Fig. 16 and Table 7. The majority (70/85) of the ex- 
tended sources are steep spectrum, with the spectral index distri- 
bution peaking at a « 0.8. The distribution of the unclassified 
sources is much broader and is less peaked, with the distribution 
stretching between a » -0.2 to a « +1. The majority of the 
sources with no information display a flat spectrum, as these are 
the sources which are too faint to be detected at the lower frequen- 
cies of 610 MHz/1.4 GHz. For these sources the value of a plotted 
is an upper limit, which explains the large peak at a ~ 0. Fig. 17 
shows Rgmrt against al l', for all sources matched to the GMRT 
catalogue. It is evident from this plot that, as expected, a greater 
number of the extended sources have a steep spectrum. 

The distribution for extended sources displays a flat spectrum 
tail, 15 sources having a£ s J < 0.5. It is likely that these extended, 
flat spectrum sources are extended sources with a flat spectrum core 
which dominates at high frequency. We might therefore expect that 
the lower frequency images of most of these sources would dis- 
play a dominant core surrounded by some, possibly fainter, exten- 
ded emission, whereas the steep spectrum extended sources would 
have relatively more pronounced lobes. Examination of the images 
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Figure 16. The spectral index distribution for sources classified as extended 
and unclassified, and those with no information. The three histograms are 
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which explains the apparent peak at a « 0. 



10' 



10" 



S/N > 20 
S/N < 20 
- R = 1 .35 
R = 2.00 



++ ++ 



++>+ . 



-0.5 



U 0.61 



0.5 



1.5 



Figure 17. Rgmrt as a function of al ^ for all sources in the GMRT area. 
Sources with peak flux / noise < 20 in the GMRT catalogue are shown 
separately and the horizontal lines at R = 1.35 and 2.00 indicate the cutoffs 
in R used when classifying the extended sources, see text for details. 



shows this to be the case for most sources; typical examples of steep 
and flat spectrum extended sources are shown in Fig. 18. 



4.4 Correlations between size and spectral index 

In Section 4.3 it is shown that the majority of the flat spectrum 
sources are not classified as extended and therefore probably com- 
pact. However, a small but significant number (15 out of 163) of 
the flat spectrum sources are clearly extended, indicating a variety 
of source types. There are also a significant number of sources (34) 
for which we have no information about their structure. 

Recent work by Massardi et al. (201 la) using the AT20G sur- 
vey found that for a? there was a clear divide between the exten- 
ded, steep spectrum sources and the compact, flat spectrum sources. 
However, when looking at the higher frequency spectral index, af, 
they found that the extended sources displayed a much broader 
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Figure 18. GMRT images of some example extended sources. The contours are drawn at (±2 V2™, n = 0, 1...7) X x mJy where x = 0.16 for (a), 0.15 for (b), 
0.086 for (c) and 0.076 for (d). The X marks the position of the 10C source and the + mark the positions of the components listed in the GMRT catalogue for 
each source. Sources (a) and (b) are steep spectrum sources with arj^ = 0.93 and 1.08 respectively. Sources (c) and (d) have flatter spectra with ajg] = 0.49 
and 0.34 respectively. 



range in spectral index. This is consistent with the flat spectrum 
tail in the spectral index distribution for extended sources observed 
here and supports the idea that a flat spectrum core is becoming 
increasingly dominant at higher frequencies. A study by Prandoni 
et al. (2006) investigated the properties of a sample of 1 1 1 sources 
using data at 1 .4 and 5 GHz and found that nearly all the extended 
or multiple sources were steep spectrum. The lack of flat spectrum 
extended sources in the Prandoni et al. study is probably due to the 
lower frequencies used, as at these frequencies the flat spectrum 
core has not yet become dominant. 



4.5 Effect on the source counts 

The 15.7-GHz source count derived from the full 10C survey is 
presented in AMI Consortium: Davies et al. (2011). The function 
fitted to the source count is a broken power law, as shown in equa- 
tion 1, 

_ dN ( 24(f)" 2 ' 27 Jy" 1 sr 1 for 2.8 < S < 25 mJy, 
dS I 376 (^) Jy" 1 sr 1 for 0.5 < 5 < 2.8 mJy. 

It is significant that the break in this power law occurs at 
2.8 mJy i.e. at approximately the same flux density as the change 
in the spectral index distribution observed here. 

To examine this further, the source counts for the complete 
sample of sources (made up of samples A and B) studied in this 
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Table 8. Source counts for the complete IOC sample from the Lockman 
Hole fields (made up of samples A and B). Steep and flat spectrum sources 
(with ffigj > 0.5 and <*Jgj < 0.5 respectively) are shown separately. Each 
flux density bin corresponds to two of the flux density bins used in present- 
ing the full IOC source counts in AMI Consortium: Davies et al. (2011) 
(except for the highest flux density bin). 



Bin start 


Bin end 


No. of sources 


No. of sources 


Area 


/mJy 


/mJy 


with a > 0.5 


with a < 0.5 


/deg 2 


9.000 


25.000 


4 


4 


4.64 


2.900 


9.000 


22 


11 


4.64 


1.500 


2.900 


27 


10 


4.64 


1.000 


1.500 


18 


24 


4.64 


0.775 


1.000 


5 


11 


1.73 


0.600 


0.775 


5 


11 


1.73 


0.500 


0.600 


6 


12 


1.73 



paper are presented in Table 8. The counts for steep and flat spec- 
trum sources are shown separately; the flux density bins are double 
the bins used in presenting the full 10C counts in AMI Consortium: 
Davies et al. (201 1). No attempt is made to fit the source count due 
to the small number of sources in each bin. However, it is clear that 
the source counts are significantly different for the two populations, 
with, as expected, a greater proportion of flat spectrum sources in 
the fainter flux density bins (5 i 5G hz ;£ 1 mJy), while the steep spec- 
trum sources dominate at the higher flux densities. Taken with the 
flattening of the overall counts for flux densities < 2.8 mJy, this in- 
dicates that the source counts for the steep spectrum sources must 
be flattening significantly at flux densities below about 1 mJy. 

There is evidence from other studies that the FRII population 
is dropping out at around the flux density where we observe the 
flattening of the steep spectrum source count. For example, Gendre 
& Wall (2008) used a sample of Combined NVSS-FIRST Galax- 
ies (CoNFIG) to construct a 1 .4-GHz source count for FRI and FRII 
sources. They found that the FRII population is dropping out as flux 
density decreases below approximately 20 mJy at 1.4 GHz, leaving 
a source population dominated by FRI sources at lower flux densit- 
ies. This change in the population at 5 i.4ghz ~ 20 mJy corresponds 
to 5 ~ 3 mJy at 15 GHz for a steep spectrum source (a = 0.75). The 
simulated source counts produced by Wilman et al. (2008) (which 
are discussed in more detail in Section 6) also show that the FRII 
population drops out at a few mJy at 18 GHz. It therefore seems 
likely that the changes in the spectral index distributions at around 
1 mJy in our 15 GHz sample are due in part to the disappearance 
of the FRII sources. However, as discussed in Section 6, the dom- 
inance of a flat spectrum population at flux densities below 1 mJy 
is in clear disagreement with the models of Wilman et al., de Zotti 
et al. and Tucci et al.. 



5 COMPARISON WITH SAMPLES SELECTED AT 
1.4 GHz 

The source population at 1 .4 GHz has been much more widely stud- 
ied than the higher frequency population and models of the faint 
population at higher frequencies are often extrapolated from this 
lower frequency data. It is therefore useful to see how the spectral 
index distribution of the 10C source population compares to that 
for sources selected at 1 .4 GHz. 

Two samples of sources at 1.4 GHz were used; the first sample 
was selected from the FIRST catalogue (sample P) and the second 
from the NVSS catalogue (sample Q). The FIRST catalogue is 




0.5 1 1.5 

log(1 .4-GHz flux density / mJy) 



2.5 



Figure 19. The flux density distribution of the sources in sample P, selected 
from the FIRST catalogue. Sources with and without a match to the 10C 
catalogue are shown separately. 



deeper so sample P provides more information about the faint 
source population; the NVSS survey on the other hand has a beam 
size comparable to the 10C survey so provides more reliable spec- 
tral indices. 

Some consideration needs to be given to the limiting flux dens- 
ities of the FIRST and NVSS catalogues and the 10C catalogue 
which they are being matched. A source with a 1. 4-GHz flux dens- 
ity of 3.4 mJy (the completeness limit of NVSS) should be detected 
in 10C unless it has a spectral index greater than 0.8. It is therefore 
expected that counterparts will be found at 15.7 GHz for the major- 
ity of the NVSS sources. FIRST, however, has a completeness limit 
of 1 mJy; a source with a 1. 4-GHz flux density of 1 mJy could be 
detected at 15.7 GHz if it had ct| 5 4 7 < 0.3. We therefore expect a 
larger proportion of lower limits on the spectral indices for FIRST 
sources. 



5.1 Sample selected from FIRST 

All FIRST sources in the 10C deep fields (the areas containing 
Sample A) were selected (see Fig. 1) giving a sample of 127 sources 
selected at 1.4 GHz. These sources were matched to the full 10C 
catalogue using a match radius of 15 arcsec, as described in Sec- 
tion 2.2. The difference in resolution between the two catalogues 
meant that there were several sources which were resolved into 
multiple components in FIRST but unresolved in 10C. For those 
cases where several FIRST sources matched to one 10C source, the 
flux densities of the individual FIRST sources were combined, giv- 
ing a sample of 105 FIRST sources, 70 of which have a match to 
IOC. 

10C maps of the unmatched FIRST sources were examined 
by eye. In ten cases, a source was visible in the 10C image - how- 
ever its flux was clearly below the 10C completeness limit; the flux 
density of such a source was found using TVSTAT in AIPS. For 
the remaining unmatched sources, an upper limit of three times the 
local noise was placed on the flux density. This allows a lower limit 
to be placed on the spectral index. 

The flux density distribution of the 105 sources in sample P is 
shown in Fig. 19. 
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Figure 20. The flux density distribution of the sources in sample Q, selected 
from the NVSS catalogue. Sources with and without a match to the 10C 
catalogue are shown separately. 



5.2 Sample selected from NVSS 

In order to have enough NVSS sources to be able to draw statist- 
ically significant conclusions, sources were selected from the deep 
areas of four 10C fields. Two of these are the fields in the Lock- 
man Hole studied in this paper, the other two fields are centred 
on 17 h 33 m , +41°48 m and 00 h 24 m , +31°52 m . There are 292 sources 
in these four fields, giving a sample of comparable size to the 
sample of 10C sources studied in this paper. The NVSS sources 
were matched to the 10C catalogue, with limits placed on the un- 
matched sources as described in Section 5.1. This sample (sample 
Q) contains 292 NVSS sources, 223 of which have a match to the 
10C catalogue. The flux density distribution of sample Q is shown 
in Fig. 20. 



5.3 Spectral index distribution of samples selected at 1.4 GHz 

The spectral index ct[ 5 4 7 was calculated for all sources in samples 
P and Q with a match in IOC. For the unmatched sources, a lower 
limit was placed on the spectral index using the upper limit from 
the 10C map. Fig. 21 shows a comparison of the spectral index 
distributions of the sources selected at 15.7 GHz and those selec- 
ted at 1.4 GHz. As expected the distributions are noticeably differ- 
ent; the sources selected at 1.4 GHz show one peak at ff} 5 4 7 as 0.7, 
while the sample selected at 15.7 GHz displays an additional peak 
at aj 5 4 7 as 0.3. The additional population of flat spectrum sources, as 
expected, are poorly represented by selecting at 1.4 GHz. This is the 
challenge that extrapolating from lower frequencies to predict the 
high frequency radio population presents. It relies on accurate mod- 
elling of the source population and of how the spectral behaviour 
of sources varies with frequency; it is particularly unreliable when 
there are no observations at nearby frequencies at the required flux 
density. The 10C sample is compared to one such model in Section 
6. 



6 COMPARISON WITH THE SKADS SIMULATED SKY 

Wilman et al. (2008, 2010) produced a semi-empirical simulation 
of the extragalactic radio continuum sky which contains as 320 rnil- 
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Figure 21. The spectral index distributions for sources selected at 1.4 GHz 
from NVSS (sample Q, shown in the top panel) and FIRST (sample P, 
middle panel) and IOC sources selected at 15.7 GHz (bottom panel). Lim- 
iting spectral indices for the unmatched sources are shown in white - note 
that these are lower limits in the FIRST and NVSS cases, therefore could 
only move to the right, and upper limits in the IOC case. 



lion sources. This simulation covers a sky area of 20 x 20 deg 2 out 
to a cosmological redshift of z = 20 and down to flux density limits 
of 10 nly at 151, 610 MHz 1.4, 4.86 and 18 GHz. The sources in 
the simulation are split into six distinct source types: radio-quiet 
AGN (RQQ), radio loud AGN of the Fanaroff-Riley type I (FRI) 
and type II (FRII), GHz-peaked spectrum (GPS) sources, quiescent 
starforming and starbursting galaxies. These simulated sources are 
drawn from the observed or extrapolated luminosity functions. In 
order to produce a sub-sample of simulated sources comparable 
to the sources observed in this work, sources with a flux density 
greater than 0.5 mJy at 18 GHz were selected from the simulation. 
This produced a simulated sub-sample of 16235 sources (sample 
S 3 ). The spectral index between 610 MHz and 18 GHz was calcu- 
lated from the flux densities in the catalogue. 

The S 3 sample is dominated by FRI sources, which make up 
71 percent of the source population (Table 9) in agreement with 
the discussion in Section 4.5. The second largest source population 
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Figure 22. The spectral index distribution of the different source population in the S 3 sample with S ig GHz > 0.5 mjy. The left panel shows FRI and FRII 
sources, the middle panel shows radio quiet AGNS and GHz-peaked sources and the right panel shows starburst and quiescent starforming soures. 
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Figure 23. Comparison of the 10C and the simulated samples. The spectral index distribution (left) and flux density distribution (right) for sources in the S 3 
sample with Sis GHz > 0.5 mjy and for 10C sources with S157 GHz > 0.5 mjy is shown. Note that for the S 3 sample the values plotted are urjgj and 18-GHz 
flux density while for the 10C sources a' 5 : 1 , and 15.7-GHz flux density are plotted. We do not expect this to make any difference to the results. Error bars show 



the poisson errors. 
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Figure 24. Comparison of the spectral index distribution and flux density distribution for sources in the S 3 sample with Sis GHz > 0.5 mjy and for 10C sources 
with S15.7 GHz > 0.5 mjy and a spectral index (oi?^) § reater man 0.3 only. Note that for the S 3 sample the values plotted are urigj and 18-GHz flux density 
while for the 10C sources advj and 15.7-GHz flux density are plotted. We do not expect this to make any difference to the results. Error bars show the poisson 
eiTors. 
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Table 9. The proportion of different source types in the sub-sample of sim- 
ulated sources selected to be directly comparable to the sources observed 
in this study. This sub-sample contains 16235 sources with a flux density at 
18 GHz greater than 0.5 mJy. The range of percentages calculated from five 
regions of equivalent size to the 10C deep survey area is also given. 



Source type 


Percentage 


Range 


FRI 


71 


57-85 


FRII 


13 


6-30 


Radio quiet AGN 


3 


0-6 


GPS 


3 


0-8 


Starburst 


4 


0-11 


Quiescent starforming 


3 


1-11 



are FRII sources, while radio quiet AGN, GPS sources and star- 
forming galaxies each make only small contributions to the simu- 
lated source population. The S 3 sample covers a region of area 400 
deg 2 while the deep 10C regions from which sample A is drawn 
cover 1.73 deg 2 . In order to investigate the possible impact of clus- 
tering on these results, five regions of the S 3 sample with areas of 
1 .73 deg 2 were selected at random. The fraction of different sources 
types in each of these five regions was calculated and the range of 
values this gave is shown in Table 9. Although the exact percent- 
ages of source vary across the five regions, the general proportions 
remain the same, indicating that clustering will not have a major 
effect on the 10C results. 

The spectral index distributions of the different groups of 
sources in the simulation are shown in Fig. 22. All FRI, FRII and 
GPS sources have been modelled assuming their extended emission 
has a constant spectral index of 0.75, hence the prominent peaks at 
a = 0.75. An orientation-dependent relativistic beaming model is 
used to find the contribution of the flat-spectrum core to the over- 
all emission from each source; this gives the flatter spectrum tail 
in each spectral index distribution. The radio-quiet AGN have been 
assumed to have a constant spectral index of 0.7. The spectra of 
the starburst and starforming galaxies have been modelled using 
thermal and non-thermal components and in the case of starbursts 
a thermal dust component has been included. 

Fig. 23 shows a comparison of the spectral index and flux 
density distributions of the observed 10C sources (sample A, com- 
plete to 0.5 mJy) and S 3 sample. It is clear that the simulation fails 
to reproduce the spectral index distribution of the 10C sources. The 
discrepancy for a > 0.7 is due to the input assumption of the model 
that all sources have a = 0.75; however the model fails to repro- 
duce the distribution for a < 0.7 with a conspicuous absence of 
sources with a < 0.3. There are two main possibilities here; either 
the distribution of sources has not been modelled correctly and FRI 
sources do not dominate at this frequency and flux density level, but 
instead a new population with flat spectra is becoming important, or 
the emission from FRI sources has not been modelled correctly and 
that their flat spectrum cores are more dominant than predicted by 
the model in this frequency range. It is also possible that starburst 
galaxies may be causing this flattening in spectral index, although 
this is unlikely as it would require the contribution of starbursts to 
be greater than modelled by at least a factor of ten. 

The actual and simulated flux density distributions are similar, 
but the 10C distribution contains a larger proportion of sources with 
flux densities less than 1 mJy. To test the possibility that there is a 
population of faint flat spectrum sources observed here which are 
missing from the simulation, the spectral index and flux density 
distributions were replotted, this time excluding all sources in the 



10C sample with a < 0.3 (Fig. 24) as essentially no sources with 
a < 0.3 are predicted by the model. The distributions are now more 
similar. 

This analysis indicates that the extrapolations of the luminos- 
ity functions coupled with the models for the effects of beaming 
on the spectra of the radio-loud AGN used in this simulation have 
failed to reproduce the observed properties of the high frequency 
population. It is worth noting that the recent models of the source 
population at 15 GHz by de Zotti et al. (2005) and Tucci et al. 
(201 1) also fail to reproduce the observed source count at flux dens- 
ities < 10 mJy, significantly under- predicting the observed number 
of sources. The updated version of the model of the 15-GHz source 
count by de Zotti et al. (2005), extracted from their website 3 , shows 
that steep spectrum sources outnumber flat spectrum sources until 
the flux density drops below approximately 2 ujy, in clear disagree- 
ment with the counts in Table 8. The de Zotti et al. model predicts 
that at 5 isghz = 1 mJy steep spectrum sources outnumber flat spec- 
trum sources by nearly a factor of three. The results in Table 8 
show that for 5 i 5G hz < 1 rnJy there are twice as many flat spec- 
trum sources as steep spectrum sources. The recent high frequency 
predictions of the source counts by Tucci et al. (2011) significantly 
under predict the number of sources observed at 15 GHz below 
approximately 5 mJy. This underprediction of the total number of 
sources could be explained by there being a greater number of flat 
spectrum sources at faint flux densities than are included in the 
model. These results highlight the difficulties inherent in predict- 
ing the behaviour of the high frequency radio source population by 
extrapolating from lower frequencies. 



7 CONCLUSIONS 

The radio spectral properties of 296 sources detected as part of the 
10C survey at 15.7 GHz in the Lockman Hole are investigated in 
detail using a number of radio surveys, in particular a deep GMRT 
image at 6 10 MHz and a WSRT image at 1 .4 GHz. Matches at other 
radio frequencies were found for 266 out of 296 sources, allowing 
their radio spectra to be investigated. For the 30 sources which were 
only detected at 15.7 GHz, upper limits are placed on their spectral 
indices. 

There is a clear change in spectral index with flux density - 
the median a [ 5 6 ] = 0.75 for flux densities greater than 1.5 mJy 
while the median a 1 ^ = 0.08 for flux densities less than 0.8 mJy. 
This demonstrates that there is a population of flat spectrum sources 
emerging below 1 mJy. This result is consistent with results from 
other studies of the spectral indices of sources at lower frequencies. 

The 10C source population was compared to two samples se- 
lected at 1.4 GHz from FIRST and NVSS at a comparable flux 
density. The spectral index distribution of these two samples is 
significantly different from that of the 10C sample selected at 
15.7 GHz, the flat spectrum population present at 15 GHz being 
poorly represented in the 1.4-GHz samples. This demonstrates the 
well-known problem with extrapolating from lower frequencies to 
predict the properties of the high-frequency population. 

The 10C sample was compared to a comparable sample se- 
lected from the SKADS Simulated Sky constructed by Wilman et 
al. (2008). The spectral index distributions of the two samples dif- 
fer significantly; there are essentially no sources in the simulated 
sample with a < 0.3 while 40 percent of the 10C sample have 

3 http://web.oapd.inaf.it/rstools/srccnt_tables 
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a 06i < 0.3. There is also a larger proportion of sources with flux 
densities below 1 mJy in the 10C sample than in the simulated 
sample - 57 percent of the 10C sources have a flux density below 
1 mJy compared to 40 percent of simulated sources. This indicates 
that the simulation does not accurately reproduce the observed pop- 
ulation at 15.7 GHz. We conclude that either there is a population 
of faint, flat spectrum sources which are missing from the simula- 
tion or the high-frequency radio emission of a known population 
is not modelled correctly in the simulation. If the relative contribu- 
tions of the different populations are modelled correctly, it is likely 
that the observed flat spectrum population is due to cores of FRI 
sources being much more dominant than the model suggests. 

Our unique, faint 15 GHz samples are of great value when 
investigating the faint, high-frequency source population. We will 
use optical and infrared data in the Lockman Hole to investigate 
the nature of this population of sources in a later paper. We are also 
extending this study to fainter flux densities at 15 GHz. 
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